Edge logarithmic corrections probed by impurity NMR 
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Semi-infinite quantum spin chains display spin autocorrelations near the boundary with power-law 
exponents that are given by boundary conformal field theories. We show that NMR measurements 
on spinless impurities that break a quantum spin chain lead to a spin-lattice relaxation rate 1/T^ ^^ 
that has a temperature dependence which is a direct probe of the anomalous boundary exponents. 
For the antiferromagnetic S=l/2 spin chain, we show that 1/T^ ^'^ oc Tlog^T instead of log^' ^ T 
for a bulk measurement. We show that, in the case of a one-dimensional conductor described by a 
Luttinger liquid, a similar measurement leads to a relaxation rate l/Tf ^^ oc T, independent of the 
anomalous exponent Kp. 
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Impurity effects in quantum spin chains are a subject 
of much interest since they are a sensitive probe of the 
peculiar correlations that develop in these magnetic sys- 
tems. For example it has been shown recently M that 
an external uniform magnetic field creates a staggered 
magnetization near impurities in the antiferromagnetic 
(AF) spin chain compound Sr2Cu03, in agreement with 
theoretical predictions 0. In S=l AF spin systems with 
a Haldane gap in the bulk there are effective free spins 
S=l/2 at the end of the chain that are revealed by 
the coupling to an impurity. 

One-dimensional quantum systems like the S=l/2 AF 
chain or the electron gas - the "Luttinger liquid" - 
are critical systems with an effective low-energy physics 
which is a conformal field theory. In both cases it is 
the massless Gaussian model that is at the heart of the 
correct description. When boundaries are introduced in 
such problems, then it has been shown by Cardy that 
there are boundary operators with nontrivial scaling di- 
mensions [0. This means that if we measure the spin 
autocorrelation in time (S(t) • S(0)) for the spin at the 
end of a chain, it decays with a power law which is dif- 
ferent from that of the bulk. Measurements with a local 
probe like the l/Ti relaxation rate in NMR experiments 
are sensitive to this power law. 

Conventional NMR is performed on nuclear spins that 
are coupled by hyperfine interactions to the atomic spin 
of the same atom : for example in the S=l/2 antifer- 
romagnetic spin chain Sr2Cu03 one may use ^■^Cu reso- 
nance H. While this is a local measurement, it involves 
the time autocorrelation of a bulk spin. We show in this 
Letter that by NMR measurements on nuclei belonging 
to spinless impurities that break the chain (or the ID 
conductor) one measures the anomalous boundary expo- 
nents of the boundary conformal field theory. We call 
this "edge NMR". In the case of the spin-1/2 chain we 
show that the bulk marginal operator manifests itself by 
a logarithmic factor in 1/Ti which is different from the 
bulk case. 

In the case of a spin chain, our main idea is illustrated 



in fig. 1. Impurities that are spinless but do have a nu- 
clear spin are introduced at a low concentration level in 
the spin chain. The nuclear spin of the impurity in fig. 1 
has substantial coupling only with its two neighboring 
spins if we are dealing with a magnetic insulator. The 
coupling may be dipolar or transferred hyperfine : its 
precise nature does not matter for our reasoning. When 
impurities are very dilute, the nuclear spin is surrounded 
by two essentially semi-infinite spin chains. So measure- 
ment of the relaxation rate Ti of the impurity nuclear 
spin involve only the fluctuations of the two end spins 
and these spins give the same contribution. If there is no 
exchange through the impurity, then we have exactly a 
probe of a boundary critical phenomena since the relax- 
ation rate is given by the following formula |^] : 



/ + 00 
dt 



(So(i) • So{0)h 



(1) 



In this equation A is the (hyperfine) coupling, Sq is the 
end spin of the semi- infinite chain and (. . .)t stands for 
thermal averaging. An example would be a S=l/2 spin 
chain realized by Cu^"*" ions and the impurity may be 
Zn which is a possible substitution for Cu in the cuprate 
family of materials. So Zn doping effectively cuts the 
chains. The isotope ^^Zn has a nuclear spin so NMR 
may be performed on this nucleus selectively. 

In the case of the S=l/2 AF chain we show that 
{l/TiY'^^'^ ex Tlog T is the universal low-temperature 
behavior instead of log^ " T as predicted |^,0 and mea- 
sured in the case of bulk NMR g . In the case of a one- 
dimensional conductor, one can perform the same kind 
of measurement by using an impurity which is insulat- 
ing, i.e. which breaks the chain. Similarly one has now a 
relaxation rate that scales as (1/Ti)'^'*^'^ ex T while in the 
bulk it is given by § (1/Ti)''"''' ex T -f T^" where Kp is 
the anomalous exponent that characterizes the Luttinger 
liquid behavior of the system. 

We first discuss the S=l/2 AF spin chain described by 
the Heisenberg hamiltonian : 



Impurity with nuclear spin only 

FIG. 1. An NMR l/Ti measurement performed on an im- 
purity with a nuclear spin but no atomic spin (white dot) is a 
probe of the time autocorrelation of the end spins of the spin 
chain (black dots). 
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(2) 



In this equation, we have introduced the exchange 
anisotropy A which is convenient for the bosonization 
description of the low-energy physics. In most materi- 
als that are known to be reasonably well described by a 
S=l/2 AF chain like KCuFg or SraCuOg, in fact there 
is isotropy, i. e. A = l. IfA = Owe have the XY 
chain which maps onto a system of free spinless ferniions 
in a band of dispersion E{k) = —Jcosk and two Fermi 
points at kp = ±7r/2. There are particle-hole excitations 
with small momentum transfer, i.e. staying on the same 
side of the Fermi "surface" and also excitations across the 
Fermi surface with momentum transfer 2kF = tt. In the 
spin language, they correspond respectively to ferro and 
antiferromagnetic fluctuations. In the bosonic language 
Pi we have : 



1 (~Y 



cos(v47r(/)(a;)). 



(3) 



The effective theory in the infrared limit for the boson is 
a free theory given by the following Lagrangian density : 



^0 = ^ [{dt^f - (d.,^ 



(4) 



For anisotropics A < 1, the chain is in a massless phase 
which is still described by Eq. (^ after a rescaling of the 
boson (j) -^ (t)/VK where K — K{A) is a parameter that 
contains the effects of interactions induced by z-exchange. 
At the isotropic point A = 1 we have K = 1/2. 

This massless effective theory leads to power laws in 
the spin fluctuations. The local correlations of a spin S 
in an infinite chain can be computed at nonzero temper- 
ature T with the result : 



{S^{t)S^O)) (X 



T' 



+ (- 



J.2 



sinh TiTt sinh nTt 



K 



(5) 
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T' 



sinh irTt 



f/^K ^ ( 



T 



sinh TiTt 



2K+1/2K 



(6) 



(omitting constant prefactors). In the isotropic case both 
spin correlations are equal to T/ sinh nTt and this be- 
havior leads to a relaxation l/Ti calculated from Eq. dl^) 
which is oc T°. This is essentially the behavior observed 
by Takigawa et al. m. There is an additional logarith- 
mic correction log^ " T due to the presence of Umklapp 
processes for the fermions |6). 

If we now consider a spin at the end of an open chain, 
we must supply our bosonic description with a boundary 
condition. It is a Dirichlet condition which is appropri- 
ate |11|. The condition (S^{x = 0)) = in the bosonic 
theory leads to 0(0) = v^/2\/2 by use of Eq. (|). This 
boundary condition can then be used to relate the chiral 
components 4>l and (fiR. Indeed (pR may be regarded as 
the analytic continuation of (Jjl on the negative axis by 
(j)ji{x,t) = —(/)£{— x,t) + •y7r/2-\/2. The Green's function 
can then be expressed in terms of the L field only and 
this field does not feel the existence of the boundary. 

The edge autocorrelations have now a different behav- 
ior from the bulk. We find that the uniform and stag- 
gered part of the correlation functions have the same be- 
havior at the edge, given by : 



{s^it)s^m ex 



{S+{t)S-{Q))^{ 



rj.2 



sinh^ T:Tt 



T 



sinh "nTt 



l/K 



(7) 



(8) 



where we have kept the leading behavior. This is due 
to the fact that spin operators have boundary scaling 
dimensions 1 instead of 1/2 in the bulk. In the isotropic 
case these correlations leads to contributions to 1/Ti that 
scale as T instead of T° in the bulk. So the boundary 
critical behavior influences directly the NMR relaxation 
rate. The leading behavior is thus {l/Tif^Be ^ x. 

We now turn to the influence of the marginal opera- 
tor. Indeed it is known that there are logarithmic cor- 
rections to the power-law behavior of the spin-1/2 chain 
p^ , p^ . In fact the effective theory given by Eq. (^ con- 
tains an additional marginal operator that is allowed in 
the isotropic case : it is due to Umklapp scattering of 
the fermions and can be written as cos(vl67r0). The 
complete effective theory is given by : 



>C0 = i [{dtd^f 



{d^(i>f] +gcos{y 



(9) 



The set of the two coupling constants K and g flows 
under renormalization group according to the Kosterlitz- 
Thouless equations. There is a single line that reaches 
the fixed point K* = 1/2 and which corresponds to an 
isotropic spin chain. Along this line, the flow of the cou- 
plings is marginal. This operator leads to a multiplicative 
logarithmic correction to the bulk magnetic susceptibil- 
ity, as predicted theoretically llj] and observed experi- 
mentally |15[. In the relaxation rate it leads to a correc- 
tion PI ex {logTy/^ which is a possible explanation for 



the experimental upturn of 1/Ti at low temperature [gj. 
Many numerical studies have been devoted to demon- 
strate ab initio its existence in the spin-1/2 chain. In 
fact it is responsible for the extremely slow approach to 
the thermodynamic limit p^ , p7[ |. 

This bulk marginal operator still exists in the theory 
with a boundary. In addition, some new operators are 
permitted at the edge. Note that the surface term dx4>L, 
while forbidden in the bulk of the chain, is a perfectly 
allowed operator in the boundary theory. If we now 
classify all possible marginal operators at the boundary 
fixed point, there are only two such operators dx4>L and 
exp(i\/87r(/)L). They have boundary scale dimension 1 
and thus are marginal on the edge. 

But it is easy to see that they are forbidden by the sym- 
metry of the problem. The rotational invariance around 
z-axis forbids exp(j\/87r0i) whereas it is rotation around 
X-axis that forbids dx(f>L- As a consequence, there are 
no extra marginal operators. Thus the Green's function 
defined by : 



G{t) 



edge 



(5o+(i)%(0)), 



(10) 



obeys the following edge renormalization group equation : 



ijr + Pbulk[g)-K- + ^ledgeig) 



G{t,g) 



edge 



0. (11) 



In this equation, g is the coupling of the marginal op- 
erator. Since it is a bulk operator its /3-function is that 
of the bulk : P(g) = (3buik{g)- The anomalous dimension 
ledge (ff) howcvcr is related to the short-distance singular- 
ity of the bulk marginal operator with an operator that 
lies at the boundary. Hence it is modified with respect 
to its bulk value. 

We compute jedgeig) by a Coulomb gas technique sim- 
ilar to that used by Singh, Fisher and Shankar |1^, ex- 
pressing all the Coulomb gas quantities in terms of 0l 
only. We find : 



Pbuik{g) = -27:g^ + Oig^), 



-/edge{g) = l-27rg + 0{g^). 



(12) 



(13) 



The result for 'yedge is different from the bulk result 
7(5) = 1/2 — ng/2 : the constant term is due to the edge 
scaling dimension of the spin operator which is twice the 
bulk value. The 0{g) coefficient leads directly to the 
logarithmic correction. The solution of Eq. (O) by the 
method of characteristics leads to an asymptotic behav- 
ior which includes a logarithmic correction : 



G{t) 



edge 



ex ^(logt)2. 



(14) 



This is the peculiar law of time decay for the autocor- 
relation of the end spin of the chain (in the bulk it is 



log"^'^i/i). This result holds for the staggered as well as 
the uniform part of the correlation function. We have 
checked that Eqs. (|l2|,^ are consistent with isotropy : 
(5(f(i)5'o-(0)) and (5g(t)5g(0)) have the same behav- 
ior. To treat the effect of a nonzero temperature we 
impose periodic conditions in Euclidean time with pe- 
riod 1/T. The system lives then on a half-cylinder. The 
Callan-Symanzik equation (O) is unchanged, though the 
Green's function depends upon an additional argument 
which is the temporal extent l/T. A standard finite-size 
scaling argument leads then to a solution of the form : 



G{t,Ty'^se ^_ 



rp2 



smh\-KTt) 



{\ogtf 



(15) 



Such a behavior leads to a relaxation rate {1/TiY'^^'^ oc 
T log T strikingly different from the bulk. Experimental 
check of this result seems feasible with the possible com- 
plication that if the nearest-neighbor coupling needed in 
edge NMR is dipolar then it may be that there are par- 
asitic effects due to the other neighboring chains. We 
note that, while we find the well-known doubling of sur- 
face exponents for the leading power law, this doubling 
does not extend to the power of the logarithm. 

The Luttinger liquid. We now apply the same line of ar- 
guments to the one-dimensional Luttinger liquid. There 
are many examples of physical systems including quan- 
tum wires, carbon nanotubes or one dimensional organic 
conductors that belong to the universality class of the 
"Luttinger liquid" . In such systems, the spin and charge 
degrees of freedom are decoupled in the low-energy long- 
wavelength limit and can be described in a bosonization 
framework by two bosonic fields, ipa for the spin and (j)p 
for the charge. For isotropic systems, the effective theory 
for 0CT is a free field theory like Eq. (Q). The theory for 
(j)p is also a free theory but with a nontrivial constant Kp 
which encapsulates all low-energy long-distance proper- 
ties of the system. The rate l/Ti is a probe of the spin 
fluctuations through the spin density wave (SDW) corre- 
lations {SDWit) SDW^{Q)) where the SDW operator 
is given by cJ^^cr^aCn^ (we consider the isotropic case). In 
a translationally invariant system these correlations have 
a uniform part with 5 « and a g « ikp part. They 
both contribute to (1/Ti)^""= leading to a law T -H T^" 
which involves explicitly the coupling Kp. In the case 
of edge NMR, we use the bosonized form of the SDW 
operator : 



SDW. 



staggered 



sin yj2nKp(j)p sin ytKq 



(16) 



The condition (S*^) = on the last spin of the chain 
leads to the Dirichlet boundary conditions (?!)ct(0) = and 
(/>p(0) = ^/TT/2yj2Kp. The staggered part contributes 
then to l/Ti through the product : 



(sm - 



'2tt4)„ (t) sin V2TT(t)a (0)) 0„ (o)=o 
(sin y/2^(l)p{t) sin y/2^(l)p{0)) ^ 



'0p(O) = 



vz_ 



(17) 



The spin part in Eq. (|17|) is exactly the same as in a 
quantum spin chain while the charge part is different 
due to the presence of Kp but also due to the peculiar 
boundary condition. As a consequence, when we com- 
pute these correlations each factor involves four nonzero 
contributions instead of two in the bulk (the electric neu- 
trality of Coulomb gas calculations is affected by the 
boundary). Let us compute the charge correlator be- 
tween two arbitrary points (xi,ii) and (2^2, ^2)- We 
give some details in the zero temperature case, a con- 
formal mapping can be used afterwards to obtain the 
corresponding expressions valid for low T. The left and 
right components of the charge Boson are related by 
i'pni^i't) = —4>pL{—x,t) + ^/2\/2Kp and the charge 
part Cp of the correlator Eq. (^7|) is thus equal to : 



{cos y^2'KKp{(l)pL{xi,tl) - (l)pL{-Xl,tl)) X 

X cos .j2TTKp{(t)pL{x2,t2) ~ (t'pLi~X2,t2))). (18) 

Evaluation of this correlation function leads to the fol- 
lowing expression : 



a;ia;2 1 



where I is equal to 



{h - t2)^ - JXI + X2)^ 
{tl-t2f-{xi-X2Y 



K,!2 



(19) 



(20) 



When approaching the edge of the system, i.e. x\ and x\ 
are 0{a) where a is the lattice spacing, the expression I of 
Eq. (|2y) goes to a constant and the correlator C becomes 
independent of time and the charge part in Eq. (|l7| ) does 
not contribute to the time decay. This is to be contrasted 
with the spin correlator in Eq. (^TJ) . There the boundary 
condition is simply (/)(j(0) — and thus in the expression 
corresponding to Eq. (|l9|), we find the difference I —l/I 
(with now Kp replaced by 1) that decays in time accord- 
ing to a power law ~ 1/i^ We thus obtain a law of the 
form T^/sinh irTt at nonzero temperature leading to a 
"superuniversal" behavior l/T^ ^'^ (x T independent of 
K„. 



We note that it has also been proposed recently 18-2C|] 
to use photoemission spectroscopy on the open end of a 
Luttinger liquid to probe the boundary effects. In this 
case the physical measurement is sensitive to the bound- 
ary spectral properties of the Luttinger liquid and also 
involves novel exponents that differ from those of the 
bulk. 

In this Letter, we have shown how edge NMR can be 
used as a probe of some boundary conformal field theo- 
ries. For the S=l/2 AF spin chain, the relaxation rate 
scales like T log^ T at the edge instead of T° log^^^ T. For 
a Luttinger liquid with gapless spin and charge degrees 
of freedom it scales as T. This universal behavior is in- 
dependent of the correlation exponent Kp which rules all 
ID behavior through anomalous power laws. 
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